Multi-Scaling Comparative Analysis of Time Series 
and a Discussion on "Earthquake Conversations" in California 
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Time series are characterized by complex memory and/or distribution patterns. In this letter 
we show that models obeying to different statistics may equally reproduce some pattern of a time 
series. In particular we discuss the difference between Levy-walk and fractal Gaussian intermittent 
signals and show that the adoption of complementary scaling analysis techniques may be useful 
to distinguish the two cases. Finally, we apply this methodology to the earthquake occurrences in 
California and suggest the possibility that earthquake occurrences are described by a colored (= 
'long-range correlated') Generalized Poisson model. 

PACS numbers: 91.30.Dk, 05.40.Fb, 05.45.Tp, 89.75.Fb 



Herein we introduce a method of multi-scaling com- 
parative analysis (MSCA) for the study of intermittent 
signals. We show that to distinguish between fractal 
Gaussian intermittent noise and Levy-walk intermittent 
noise the scaling results obtained using diffusion entropy 
analysis (DEA) should be compared with those obtained 
from both finite variance scaling methods (FVSM) and 
probability distribution functions (pdf) 0, 0- Finally, 
we apply MSCA to the seismic data of California and 
suggest that, instead of being described by a statistics 
according to which the waiting times between Omori's 
earthquake clusters are uncorrelated from one another, 
as the traditional Generalized Poisson model 0, 0] or a 
recent Levy- walk like model assume the data may 
also be characterized by intercluster 1/f long-range cor- 
relations that may disclose the earthquake conversations 
recently suggested by Stein jjj. 

Hurst, in his pioneering work 5], introduced the notion 
of reseated range analysis of a time series that takes the 
scaling form of (R/S)(t) oc t H (H is now called the Hurst 
exponent). This stimulated Mandelbrot to introduce the 
concept of fractional Brownian motion (fBm) 6]. In a 
random walk context the value H = 0.5 indicates un- 
correlated noise, < H < 0.5 indicates anti-persistent 
noise and 0.5 < H < 1 indicates persistent/long-range 
correlated noise |f| . Alternative scaling methods applied 
to a time series where i = 1, 2, . . . , focus on the au- 
tocorrelation function C(t) oc £i+t) oc t 2H ~ 2 , on the 
power spectrum representation Ps{f) oc / 1 ~ 2ff and 
on the evaluation of the variance of the diffusion gener- 
ated by {&} @ ( V(t) oc t 2H ). All such scaling methods 
are related to the original Hurst's analysis and yield his 
H-exponent. These techniques, referred by us as FVSM, 
assume a finite variance and according to the central limit 
theorem (CLT) the underlying statistics are Gaussian. 

Recently, Scafetta et al. introduced a complemen- 
tary scaling analysis, the DEA, that focuses on the scal- 
ing exponent S evaluated through the Shannon entropy 
S(t) of the diffusion generated by the fluctuations 



of a time series0 0. Here, the pdf of the diffusion pro- 
cess, p(x, t), is evaluated by means of the subtrajectories 
x n(t) = X) i=0 £i+n with n = 0,1,... The pdf scaling 
property for a fractal anomalous diffusion takes the form 
p(x, i) — t~ s F (x t~ S ) , and its entropy increases in time 
as S(t) = — J^gPiXjt) \n[p(x,t)] dx = A + 8 ln(i), where 
A is a constant. One can also examine the scaling proper- 
ties of the second moment for the same process using the 
FVSM. One version of FVSM is the standard deviation 
analysis (SDA) Q which is based on the evaluation of the 
standard deviation D(t) of the same variable x and pdf 
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p(x,t), and yields D(t) — y/ (x 2 ;t) — (x; t) 2 oc 

Note that the entropy S(t) does not require the vari- 
ance of the pdf p{x,t) to be finite 0- The existence of 
scaling for a process with a diverging second moment im- 
plies that DEA is complementary to and not simply an 
alternative to FVSM. So, the scaling exponent 5 is con- 
ceptually different from the Hurst exponent H measured 
by means of the FVSM. This suggests that the scaling ex- 
ponents 6 and H may fulfill multiple relations according 
to the process under study and, therefore, the combined 
use of DEA, SDA and pdf analysis may increase our un- 
derstanding of complex phenomena through a MSCA. 

Herein we focus on the statistics of intermittent noises. 
The simplest way to represent intermittent noise is 
through a dichotomous representation in which the value 
£ = 1 indicates the occurrence of an event and the value 
£ = represents no-event 0, Q, 0, 0] . An intermit- 
tent noise is characterized by the correlation properties 
of the waiting time sequence {t^} between consecutive 
events and by its waiting time distribution ip{r). There 
are two basic distinct forms of intermittent noises: 
1) Fractal Gaussian intermittent noise is character- 
ized by a long-range correlated waiting time sequence, 
(tj Ti +t ) oc t 2H ~ 2 , and by a finite variance waiting time 
distribution ?/>(r) whose form may be, for example, that 
of a Gaussian, exponential or Poisson distribution. The 
diffusion generated by a fractal Gaussian intermittent 
noise is a particular type of fBm and satisfies the asymp- 
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FIG. 1: DEA end SDA of [A] a fractal Gaussian intermittent 
noise with V , ( r ) x exp(— r/7) with 7 = 25 and H = S — 
0.75; the fractal Gaussian relation Q of equal exponents is 
fulfilled; [B] a Levy-walk intermittent noise with tf>(r) oc r _M 
and = 2.5; note the bifurcation between H — 0.75 and 
S = 0.67 caused by the Levy- walk relation @. 



totic scaling realtion between indices 
S = H . 



(1) 



We refer to as the fractal Gaussian diffusion relation. 
If the long-range correlations of {n} are destroyed via 
shuffling, the new intermittent sequence is characterized 
by the value H = S = 0.5 of random time series. Fig. 
1A shows the scaling properties of a computer generated 
fractal Gaussian intermittent noise with an exponential 
waiting time distribution and H = 5 = 0.75. 
2) Levy-walk intermittent noise is characterized by an 
uncorrelated waiting time sequence, (r, tj) oc 5ij, and a 
Levy or an inverse power law waiting time distribution 
ip(r) oc (T + t)~ m , with 2 < /i < 3 that ensures that 
although the first moment of t is finite, the second mo- 
ment diverges. The presence of a Levy-walk process in a 
given time series can be detected 0, UJ| by mean of the 
following asymptotic relation among the three exponents 



0.5 < 5= (3- 2H)~ 



rr 1 < 1 



(2) 



We refer to (5J as the Levy-walk diffusion relation. Inter- 
esting applications of this type of noise have been found 
in several Levy phe nomena including the distribution of 
solar flares 0, Till 0, 0, 0] . In the particular case 
in which the sequence {77} is correlated the scaling ex- 
ponents 6 and H are larger than the values predicted 
by Eq. Fig. IB shows the scaling properties of 

a computer generated random Levy-walk intermittent 
noise with /j, = 2.5 that has H — 0.75 and 5 = 0.67. 

We stress that the Levy-walk relation (JSJ is fulfilled 
if the waiting times {rj} are uncorrelated, in which case 
any shuffling of {n} would not alter the scaling exponents 
H and S. In fact, the super-diffusion scaling exponents 
0.5<r5<_ff<lofa Levy- walk intermittent noise are 
related to the fatness of the waiting time inverse power 
law tail, as measured by the exponent /z. Contrary to 
a fractal Gaussian intermittent noise, this Levy scaling 




does not imply a temporal correlation, or a historical 
memory, among events because the occurrence of future 
events is independent of the frequency of past events. 

We also observe that there exist particular intermit- 
tent sequences obtained by mixing Levy and Gaussian 
noises |l2j . with a Levy memory beyond memory 
by substituting an event with a cluster of events 
these cases the asymptotic properties of the scaling expo- 
nent H and 6 are expected to depend on the component, 
Gaussian or Levy, with the strongest persistence. 

The relations JTJ and (J2J, and the correlation/shuffling 
effects indicate that the DEA should be jointly used with 
the FVSM and/or pdfs. The adoption of a single tech- 
nique can lead to a misinterpretation of the character- 
istics of a phenomenon, because Levy-walk intermittent 
noise can be confused with fractal Gaussian intermittent 
noise, and uncorrelated noise of one kind of statistics can 
be mistaken for correlated noise with another kind of 
statistics. Figs. 1A and IB clearly show that the deter- 
mination of only one of the two exponents H and r5 is 
not sufficient to conclude whether a phenomenon is char- 
acterized by a Levy-walk intermittent statistics or by a 
fractal Gaussian intermittent statistics. So, we suggest a 
MSCA by combining complementary techniques. 

Recently, DEA has been applied by Mega et al. [4| to 
study the time distribution of earthquakes in Southern 
California (20°-45° N latitude and 100°-125° W longi- 
tude) from 1976 to 2002. The catalog ^ is complete 
for local events with magnitude M > 3 since 1932, for 
M > 1.8 since 1981 and for M > since 1984. The 
time intervals between large earthquakes was studied in 
Ref. Q by setting a temporal variable £(t) = 1 at the oc- 
currence of an earthquake with a magnitude larger than 
a given threshold M t , and by setting £(t) =0 when no 
earthquake of the specified magnitude occurs. We refer 
to {t{\ as the waiting time sequence between consecutive 
earthquakes with M > M t . So, the intermittent sequence 
was analyzed by means of the DEA and the measured 
scaling exponent was 5 = 0.94 ±0.01. The authors of 
Ref. |4| concluded that the time intervals, t'™|, between 
two consecutive Omori's earthquake clusters |l7j is mod- 
eled by an inverse power law ip{T^ n ~ l1 ) oc (t^)^^ with an 
exponent fi = 2.06 calculated via Eq. (J2J. This calcu- 
lation was based on the traditional assumption that 
the waiting times between such clusters are uncorrelated 

(( T l™^ r | m ') = ^*i> i m ply m g that the observed superdif- 
fusion is induced by a Levy-walk between the Omori's 
clusters. Finally, Mega et al. Q showed that a synthetic 
sequence produced with Omori's uncorrelated clusters, 
(t]™' rj m ' } = Sij , temporally distributed according to an 
inverse power law ip{T<- m >) oc (V" 1 !)^ with /1 = 2.06, 
generates a superdiffusive process with 6 — 0.94. 

However, the authors of Ref. 0] did not make the im- 
portant distinction between Levy-walk and fractal Gaus- 
sian intermittent noises such as we did above. We showed 
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FIG. 2: pdf of the waiting times Ti of earthquakes with a 
magnitude M > M t = 1, 2, 3 and 4. The initial P(r) oc 1/r 
is the Omori's law |l7|. 



that a scaling exponent in the range 0.5 < 5 < 1 can be 
associated with either a correlated fractal Gaussian inter- 
mittent noise or with an uncorrelated Levy-walk noise. 
Consequently we apply a MSCA to determine which of 
the two statistics better describes the data. 

Fig. 2 shows the waiting time pdfs between earth- 
quakes using four magnitude thresholds M t = 1, 2, 
3 and 4. The pdfs show an initial Omori's law 01 
(P(t) oc 1/r), but the pdf tails present a large inverse 
power law exponent /i > 4 and may even approach an ex- 
ponential (or Poisson) distribution asymptotically. The 
Omori's law is determined by the short-range correlated 
aftershocks 0] and lasts for a time that increases with 
the magnitude threshold. If the waiting time distribu- 
tion between Omori's clusters were an inverse power law 
with (i = 2.06 it might be expected that by increasing 
the magnitude threshold, most of the aftershocks could 
be cut off and the tail of the distribution could converge 
to an inverse power law with fi = 2.06. This does not 
seem to happen. Therefore, such an /i = 2.06 inverse 
power law, if it is real, cannot be observed in this way. 

Fig. 3 compares the DEA and SDA applied to the 
intermittent sequence of earthquakes with M > 1. 
Different magnitude thresholds give similar results. If 
these data corresponded to a random intermittent Levy- 
walk and if the curves shown in the figure corresponded to 
the asymptotic regime, the condition (J2J) interrelating the 
exponents should hold. A rigorous DEA fit in the range 
[2 7 : 2 15 ] gives S = 0.944 ± 0.008 implying a Levy-walk 
H = 0.97 ± 0.005. Instead, the same-range SDA fit gives 
H = 0.943 ± 0.004. The error analysis seems to confirm 
better the Gaussian relation of equal scaling exponents 
Q because S and H overlap within the statistical error 
as in Fig. 1A, while the difference between the measured 
H and the Levy-walk H is statistically significant (p < 
0.01). By shuffling the earthquake waiting time intervals 
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FIG. 3: DEA and SDA of the intermittent time signal £(t) 
for magnitude M > Mt = 1. The data are fitted with scaling 
exponents S = 0.944 ± 0.008 and H = 0.943 ± 0.004. The 
uppermost solid line with H = 0.97 corresponds to the ex- 
pectation of H if the Levy-walk condition holds true. 



{t.;} we get H = 0.5. Finally, by directly applying DEA 
and SDA to the waiting time series {t^}, we again get 
S = H = 0.94. These findings suggest that the data do 
not fulfill the Levy-walk relation @ and that it might 
be more likely that the Californian earthquakes are long- 
range temporal correlated according to the persistence of 
a fractal Gaussian intermittent noise with H sa 1 known 
as 1 // or pink noise Q ■ 

The curve with circles in Fig. 4 shows the DEA ap- 
plied to a synthetic earthquake catalog obtained by color- 
ing a kind of Generalized Poisson model for earthquakes. 
First we generated several Omori's clusters exactly as 
done in Ref. 0, that is, by assuming that the num- 
ber of earthquakes in a cluster follows an inverse power 
law distribution with exponent equal to 2.5 and that 
the events within the same cluster are temporally dis- 
tributed according to the Omori's law, that is, an inverse 
power law with exponent equal to p — 1. We generated 
a total number of events equal to the total number of 
earthquakes in the catalog with a magnitude threshold 
Mt = 1. However, contrary to what was done in Ref. 0], 
we do not randomly ((rj m ^ 7"j m ') cx Sij) position these 
clusters according to an inverse power law intercluster 
waiting time distribution -0(r[ m l) cx l/r** with /i = 2.06. 
Instead, we distribute the clusters according to a 1/f 
fractal Gaussian intermittent noise {i".!" 1 '} that is, with 
(t]™ 1 ' t\T}) cx t 2H ~ 2 and H « 1. The intercluster waiting 
time distribution tp(T^) cx exp(— t^/t) shown in the 
insert of Fig. 4 could be substituted with any other distri- 
bution with finite variance. Fig. 4 show that the model is 
able to reproduce the same superdiffusion pattern shown 
by the data. Finally, the curve with triangles shows the 
reduction of long-range persistency of a synthetic catalog 
obtained with the same clusters of above but temporally 
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FIG. 4: DEA applied to a 1/f long-range correlated Gen- 
eralized Poisson model (circle) for earthquakes for M > 1. 
The intercluster waiting time pdf i/;(r' m ') is exponential and 
(r| ml rj™ ] ) oc t 2H - 2 with H f* 1. The curve with triangles 
refers to the case in which the intercluster waiting time se- 
quence {t'™'} is randomized such that (i~j m ' r' m ') oc Sij. 



distributed after having shuffled, to randomize, the same 
intercluster waiting time sequence {7"j m '}. 

However, Eqs. QJ and @ are fulfilled only asymp- 
totically where the Central Limit Theorem for Gaussian 
rocesses or its generalization for Levy processes apply 
There might be the possibility that Fig. 3 as well 
as figure 2 in Ref. 4] do not show the asymptotic limit 
but a transition regime that is strongly superdiffusive 
(<5 w 1) because of the Omori's intracluster short-range 
correlations. A long transition regime is also evident in 
the curve with triangles shown in Fig. 4 that refers to 
an uncorrelated intercluster waiting time sequence and 
should asymptotically converge to 5 = 0.5. In fact, also a 
random Generalized Poisson model or an "ETAS" model 
|l8j with appropiate parameters may generate a synthetic 
catalog showing superdiffusive properties similar to the 
real data within a time range [19j. However, we observe 
that such a result may depend too strongly on the par- 
ticular parameters used in the models. 

In conclusion, we have discussed some of the difficulties 
that can be encountered in interpreting intermittent se- 
quences and shown that models with alternative statistics 
can reproduce some pattern of a time series equally well. 
This fact suggests the need of an analysis involving com- 
plementary tests. In particular we showed how to distin- 
guish fractal Gaussian intermittent noise from Levy-walk 
intermittent noise using MSCA. This methodology has 
important application in the analysis of phenomena hav- 
ing intermittent signals because different statistics imply 
different dynamics. Our analysis supports the idea that 
earthquakes generate strain diffusion, whose propagation 
over hundreds of kilometers induces remote sismic activ- 
ity 0, |2jj . This propagation according to our statistical 



analysis produces correlations in the time intervals be- 
tween earthquake clusters. In fact, the thesis that earth- 
quakes are assembled into uncorrelated Omori's clusters, 
(r] m ' rj m ' } = Sij , as both the standard Generalized Pois- 
son model and the Levy-walk model require, seems 
unrealistic. We suggest that it is more plausible that 
earthquake clusters are 1/f long-range correlated and, 
perhaps, they are subclusters of a larger Omori cluster. 
In fact, a 1/f noise can be generated by the superposi- 
tion of relaxation processes within a wide range of ener- 
gies that may well describe the coexistent stress al- 
terations caused by old and recent, as well as large and 
small shocks. Thus, the 1/f long-range correlations may 
imply that earthquake occurrences may strongly depend 
on the geological history of a vast region. 
Acknowledgment: N.S. thanks the ARO for support 
under grant DAAG5598D0002. 
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